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Abstract 



We propose a method for microscopic calculations of nuclear ground- 
state properties in the framework of density functional theory (DFT). 
We discuss how the density functional is equivalent to the effective ac- 
tion for the density, thereby establishing a constructive framework for 
DFT calculations from nuclear forces. The presented approach starts 
from non-interacting nucleons in a background potential (a simple ap- 

■ proximation for the mean field). The nuclear forces are then gradu- 
, ally turned on, while the background potential is removed. The evolu- 

^ I tion equation yields the ground-state energy and density of a system of 

■ interacting nucleons, including exchange-correlations beyond the RPA 
, approximation. The method can start from non-local low momentum 

(Mowfc) or chiral interactions. 

! 1 Introduction 

Microscopic approaches to large nuclei are inherently difficult, since nucleons 
^ ' interact strongly and the number of possible configurations increases rapidly 

' with the number of nucleons. As a consequence, for the largest nuclei the 

computational method of choice has been mean field theory. A promising 
framework to extend and systematize mean-field calculations is density func- 
tional theory (DFT). DFT is the underlying theory in modern nuclear structure 
approaches to heavy nuclei, see e.g., [1-6]. 

In this talk, we discuss an approach complementary to parameterizing the 
nuclear density functional. We present a DFT method to calculate nuclear 
ground-state properties from microscopic two-nucleon and few-nucleon forces. 
In a DFT framework, the dynamic degrees of freedom are one-body densities, 
which simplify a microscopic description for larger systems and are compu- 
tationally feasible compared to solving for the many-body wave function. In 



1 



addition, DFT can start from the softer low momentum (VJowfe) or chiral in- 
teractions [7-13], which are non-local. This is particularly important as the 
effective field theory approach offers at present the only consistent and prac- 
tical counting scheme to construct three-nucleon forces. We start by showing 
that the effective action formalism provides a microscopic foundation for DFT. 
We focus especially on how exchange-correlations are included. A direct calcu- 
lation of the effective action for non-interacting fermions shows possible expan- 
sions of the density functional when interactions are included. We then present 
a functional RG-inspired method to calculate nuclear ground-state properties 
from microscopic interactions and discuss possible extensions. 

2 Effective action formalism and DFT 

We consider non-relativistic fermions in a time-independent background po- 
tential V{'x) (e.g., the ionic potential for electronic systems; for nuclei, V — 0) 
interacting through two-body forces U (x, y). We note that three-body interac- 
tions can be easily included and for simplicity we also suppress isospin indices. 
The action S consists of a one-body Si and an interacting two-body part 5*2 
(in Euclidean space-time with spin indices a and h = 1) 

(T 

(1) 

The generating functional W[K] for connected density correlation functions 
is given by a path integral over all configurations^ 

^w[K] ^ j p[^t]p[^] e-5[^t,^]+(^t^).K^ (2) 

where we have introduced spin-dependent sources K„^„i {x) coupled to the pos- 
sible density operators 'ipl^^x) ipfji{x). The strategy for introducing the source 
term is (loosely) analogous to the calculation of the magnetization of spins on 

^For simplicity, we introduce the short-hand notation ip^ ■ ip = J2a I i'ti^) i^aix) or 
{il>^il>) ■ M = J2a a' I dxtpl{x) tpcr' {x) Ma^a' {x) ■ This allows to rewrite the action in compact 
form 5*1 [t/j'I' , V'] = "0^ ' ■ -0 and 5*2 [V"^ , V"] = 5 (V'^V') ' U ■ {tp^tjj), where the free propagator 
G includes the background potential V. 
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a lattice. For spin systems, one introduces a magnetic source H coupled to the 
total spin of the system, H Si- The variation of the source H then probes 
all different configurations and allows one to determine the spin configuration, 
i.e., the magnetization, which minimizes the free energy. For ground-state 
properties of many-body system, the source is coupled to the density operator 
and variations probe different density profiles. Thus, by varying the generating 
functional iy[-ft'] with respect to the source K one obtains the density Pa,a'{x) 
in the presence of the source (and the ground-state density for K = 0) 

Assuming that Eq. (jS)) is invertible, it can be used to express the source in terms 
of the density K[p], and one can perform a functional Legendre transformation 
from to the effective action T[p] given by 

r[p] = -W[K] +K-p. (4) 

The effective action has the advantageous property that it is minimal at the 
physical (zero source) ground-state density p = Pg 
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6T[p] 



6p 



0. (5) 



P=Pe 



Moreover, the effective action at the minimum yields the ground-state energy 
of the system Egg. This follows from the spectral representation of the path 
integral = Yin e~^^" and for K = one has in the zero temperature limit 

Egs = E[pg,] = lim ^ r[pg,]. (6) 

Summarizing, the effective action is a functional of the density which is minimal 
at the ground-state density and reproduces the energy. Therefore, the effective 
action for the density is equivalent to the energy density functional in the sense 
of Hohenberg-Kohn jT3] (more precisely for time- independent sources). This 
has been pointed out in the literature, see e.g., recently [15-18].^ One advantage 

^Two more comments are in order. First, a chemical potential corresponds to a constant 
source, which requires some attention in the calculation of the effective action and leads to 
a minimum condition ST[p]/Sp\p=p^^ — ^ for homogeneous systems, since the source and the 
chemical potential cannot be distinguished. Second, while the densities are to be expanded 
in a set of orbitals, Pgs(x) — J2n Cri|0n(x)p, density fluctuations such as p{x) — Pgsix) must 
be expanded in a set of basis functions whose integral vanishes. 
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of the effective action formalism is that it presents a constructive framework 
for microscopic calculations of the density functional. For completeness, we 
further prove the Hohenberg-Kohn universality (i.e., the trivial V dependence) 
of the density functional. This simply follows by noticing that the one-body 
potential V enters in the many-body dynamics in the same way as the external 
source K. Therefore, the dependence of the generating functional on V 

is linear, = PFy=o[-^ ~ V] and the Legendre transformation gives 

T[p] = -W[K]+K-p = -Wv=o[K-V] + {K-V + V)-p = Tv=o[p] + V-p, (7) 

where K = K — V is taken as source for the Legendre transformation. 



3 Non-interacting fermions 

After discussing the microscopic basis of DFT using effective action techniques, 
we consider the non-interacting case. This provides insight into physical ex- 
pansions of the effective action and how exchange-correlations will be included. 
In the non-interacting case, the generating functional is a Gaussian integral 
over the fermion fields and one has 



mrec [K] = Tr log[ G-'-K] = Ti log G"^ - ^ Tr [{G ■ Kf] . (8) 

n=l 

The density in the presence of the source is given by 6W{^cc[K]/SK^,j/{x), 
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p^y (x) = - G^,^^{x, x)-}_^ G^>Jx, ..)-K-{G- Kf ■ x), (9) 

n=0 

where the source is diagonal in the space-time index x and summed over indices 
are indicated by .. . For = 0, we have the familiar expression for the ground- 
state density, Pgs;o-,cr' (^) = — G a' ,cr{,x , x) . A second variation of the generating 
functional yields the particle-hole propagator G (or density-density correlator). 
From Eq. Q, we find for vanishing source in the ground state 



(JX,Y = '^{ai,a[,x),{a2,cr'2,y) 



(5^^2,4 (y) SK^^^^n 



X] 



K=0 



= {^pl^{x) ^p^'^{x) ipl^iy) ^jj^'^iy)) = - G^'^^^^{x,y) G„'^^„^{y,x), (10) 

where we have introduced a combined index X = (o"i, a[,x). It is important to 
realize that the second functional derivatives of the generating functional and 
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the effective action are related (as for tliermodynamic potentials conjugate by 
a Legendre transformation). Thus, one has 

y _ s'wjK] _ ^^^^ 

6px 5pY J 5Kx SKy 

Eq. (fTTj) relates the dressed particle-hole propagator (for the interacting sys- 
tem) to the curvature of the effective action in the ground state. We will later 
see that all exchange-correlations are built up in this way. 

For the Legendre transformation, we need the inverse expression to Eq. Q, 
i.e., the source as a function of density. By rewriting Eq. (jH)) in compact form 

Px = p,s;x + Gx,..-K„ + 0{K'), (12) 

it is clear that the source K, and thus the generating functional 14^ [i^] can be 
expanded in a power series in p — pgs around the ground-state density. For the 
source, one similarly has 

Kx = -G^l ■ (P - Pgs).. + 0{ip - Pgs)'). (13) 

We can now perform the Legendre transformation order-by-order, where 
the effective action has the same expansion around the ground-state 

oo „ ^ 

free[p] = T'^I + J2 ^t^L,Xu...,X„ iP " Pgs)^i " " " (P " Pgs,A)x„, (M) 



with ground-state energy E^s = hm r[°gg//3. Moreover, by keeping the X- 

dependence of the expansion coefficients, no local density approximation has 
been made. A straightforward calculation of the expansion coefficients yields 

r;°l = -TrlogG-\ (15) 

^iTee;X,Y = ^xV' ^^^^ 

^i2e;Xi,X2,X3 = / Syi,Y2,Y3'-^Yi,Xi^Y2,X2^Y3,X3^ (17) 

JYi,Y2,Y3 

where S^^^ denotes the totally symmetric 3-propagator ring 
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^xlx2,X3 = ■^{Ga[,a2{xi,X2)G„'^^„^{x2,X3)G„'^^„^{x3,Xi) + permutatlous) . 

(18) 

Higher-order expansion coefficients have a similar structure as T^^\ With the 
expansion around the ground state, we next present a method that includes the 
interactions among nucleons in a non-perturbative manner. It is reminiscent 
of the quasiparticle picture and the adiabatic turning on of the interaction. 
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4 Evolution equation for the density functional 



We are interested in solving for the ground-state properties of the physical 
system, where the nucleons interact by means of the microscopic interactions, 
without the presence of a background potential. To this end, we introduce a 
control parameter A, with ^ A ^ 1, by replacing 



where we allow the background potential V to depend on A, as long as the cor- 
responding eigenf unctions are localized and physically reasonable. At A = 0, 
when the two- and higher-body interactions are absent, the system is trivially 
bound by employing an attractive background potential (a simple guess for the 
mean field). The nucleons fill up the lowest-lying states to the Fermi energy. 
(Generally, the initial background potential should be chosen to improve the 
convergence of the presented algorithm.) For A = 1, we have the physical 
system, where the background potential is switched off. We use the control 
parameter to evolve a simple shell-model configuration to the fully-interacting 
system including many configurations. This is depicted in Fig. ^ We note 
that such an RG-inspired method does not require the coupling strength to be 
small. The smallness of the interaction in perturbation theory is taken over 
by the differential change, 6X serving as a "small parameter" . 



Figure 1: A graphical illustration of the initial condition of the evolution 
equation (left figure) and the final, physical system of nucleons (right figure). 
The evolution equation interpolates between the two cases. 

As the control parameter is gradually increased from A = to A = 1, the 
density functional follows an evolution equation given by [T^ 



(19) 





(20) 
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where the derivation is similar to ^7]. Since Eq. ()20j) is exact, the solution 
for general density configurations is cumbersome and is not directly related 
to physical quantities. Therefore, we propose to solve the evolution equation 
only in the vicinity of the physical point, i.e., expanded around the evolving 
ground-state density Pgs.A- Before we discuss further details of the approach, 
we separate off the background and Hartree contributions to the density func- 
tional, since their density dependence is known exactly. Thus, we introduce the 
kinetic and exchange-correlation part Tx[p] of the density functional through 

T,[p] = il-X)V,-p+^p-U-p + f,[p], (21) 

which leads to a simplified evolution equation for Tx[p] 

As the ground-state density changes under the RG, we expand the effective 
action around the current ground state pgs,A of the evolving system. This is 
the underlying theme and efficacy of our RG-inspired approach. As in the 
non-interacting case, we have 

TaIp] = rf + 5^ / - rj;!^,...,^^ ■ (p - p^,,)^, ■ ■ ■ (p - p,,,,u„. (23) 

We make the approximation and truncate the expansion for Fx at a given order 
N-p in density fiuctuations. Since one expects five- and higher-body correlations 
to be relatively small, due to the Pauli principle, a truncation at A^r = 3 or 
N-p = 4 seems a reasonable choice. Obviously, the convergence of such an 
expansion must be checked by including higher-order contributions. (We note 
that, as for the Fermi liquid RG approach used in [201121], the fiow equation 
generates upon integration an infinite set of diagrams in the underlying two- 
body interaction U.) 

The truncation of the evolution equation at three-density correlations. 
At = 3, leads to coupled RG-like equations for the expansion coefficients, 
i.e., for the kinetic and exchange-correlation (kxc) energy -Ekxc, the evolving 
ground-state density, the dressed particle-hole propag ator Gx = {T^x^ + XU) ^ 
and the three-density correlator . To provide further insight, we focus on 



6p 6p 



XU 



(22) 
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Figure 2: Graphical representation of the RG-hke equation for the ground- 
state energy -Ekxc- -E'kxc, the ground-state density pgs,A and the density- density 
correlator G\ are denoted by a blob with n = 0,1,2 dots respectively, which 
label the index structure. (1 — A)Va is represented by a wiggly line and the 
two-body interaction t/ by a dashed line. Open dots further label open indices, 
and filled dots represent indices which are traced over. 



the fiow equation for the ground-state energy and density, given by 



P dxEi,^c[pgs,x] = dxT 



(0) 
A 



-il-X)Vx-dxPgs,x->^Pgs,x-U-dxPgs,x + ^ Tr 



U-Gy 



(24) 

Here, the contributions are due to the change in the ground-state density (first 
and second terms in Eq. ((211)) and due to the dressed particle- hole propagator 
Gx (last term). The fiow equation is shown diagrammatically in Fig. |21 The 



U-G> 



Y2 



Uyi,Y2 Gx;Y2Xi- 



and 



last term in Eq. (j24|) explicitly reads Tr 

thus to lowest order is the Fock contribution to the energy. This establishes 
that all exchange-correlations are built in through the dressed particle-hole 
propagator Gx- (Obviously, such a theory with densities as variables is not 
mean-field theory.) Moreover, we note that one recovers the standard RPA 
approximation when only Gx is maintained, i.e., Nr = 2 for the effective 
action. The fiow equation for the ground-state density is given by 



dxPgs,X;X = -dxlil - A) Vx] ■ Gx-x - pgs,x ■ U ■ Gx-x 



U -Gy 



rlf ^ ■ Gx;x ) ■ Gx 



(25) 



where the expansion coefficients are totally symmetric and the indices not 
written explicitly are traced over. In the first iteration of the fiow equation, F*^^) 
is given by the 3-propagator ring with three inverse density-density correlators, 
which cancel with the three G in the last term of Eq. ()25|) . The remaining 

exchange contributions from the 3-propagaptor ring U ■ S 



(3) 
A;X 



will lead to 



a quasiparticle-like distribution over the initially unoccupied orbitals. 
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5 Conclusions and extensions 

In summary, the effective action formalism is a constructive framework for tlie 
density functional. It allows for calculations of ground-state properties starting 
from microscopic two- and three-nucleon interactions, which can be non-local 
such as the soft low momentum (V^owfc) or chiral interactions [7-13]. The main 
advantage of DFT lies in using densities as degrees of freedom, as compared to 
many-body wave functions. As a consequence, microscopic DFT approaches 
should scale well to larger nuclei. There is some insight from electronic systems, 
where DFT (see e.g., |22j) and the coupled cluster method give overall good 
results for ground-state properties in larger systems. 

We believe that the presented DFT approach may be a step towards the 
derivation of nuclear structure of heavier systems from nuclear forces. We are 
currently working on the feasibility of this approach for strongly-interacting 
one-dimensional models ^Hj- Systematic, microscopic approaches are particu- 
larly important for extrapolations to neutron/proton-rich systems, e.g., due to 
the constraints on the isospin dependence from nuclear forces or to maintain 
a complete basis of induced non-central interactions |23j . 

Extensions of this method are under investigation. A projection on vanish- 
ing center-of-mass kinetic energy can be implemented by introducing sources 
for the center-of-mass motion. The harmonic oscillator frequency of the back- 
ground potential can be adapted under the evolution to improve the conver- 
gence (reflecting the variational character of the single-particle basis). Explicit 
pairing correlations can be introduced by coupling a source to the off-diagonal 
densities ■ipa{x)'ipa' (x) (and the Hermitian conjugate). The additional minimum 
condition of the effective action is then equivalent to the BCS gap equation. 
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